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Abstract 

Double differential cross sections between 0°-12° were measured for the ^'^Zr(n,p) 
reaction at 293 MeV over a wide excitation energy range of 0~70 MeV. A multipole 
decomposition technique was applied to the present data as well as the previously 
obtained ^''Zr(p, n) data to extract the Gamow- Teller (GT) component from the 
continuum. The GT quenching factor Q was derived by using the obtained total 
GT strengths. The result \s Q = 0.88 ± 0.06, not including an overall normalization 
uncertainty in the GT unit cross section of 16%. 

Key words: Charge exchange reaction; Gamow- Teller strength; Gamow- Teller sum 
rule 
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The (p, n) reaction at intermediate energies (Tp > 100 MeV) provides a 
highly selective probe of spin-isospin excitations in nuclei due to the energy 
dependence of the isovector part of nucleon-nucleon {NN) t-matrices [1]. The 
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0° spectrum of this reaction is marked by the dominance of the Gamow- Teller 
(GT) giant resonance (GTGR), which is the aT± mode [2,3]. There exists 
a model-independent sum rule, Sp- — S0+ = 3{N — Z), where S'/3± is the 
total GT strength observed for the type [4]. Surprisingly, however, only 
a half of the GT sum rule value was identified from (p, n) measurements in 
the 1980's on targets throughout the periodic table [5]. This problem, the 
so-called quenching of the GT strength, has been one of the most interesting 
phenomena in nuclear physics because it could be related to non-nucleonic 
degrees of freedom in nuclei. At that time it was often assumed that the 
missing strength was shifted to the energy region of the A excitation due to 
coupling between nucleon particle-hole (ph) and A-isobar nucleon-hole (A/t) 
states [6,7]. However, part of the quenching is due to the nuclear configuration 
mixing between Iplh and 2p2h states [8,9]. To discuss the contribution of 
2p2h states quantitatively, one should search for missing GT strength in the 
continuum, excitation energy region of 20-50 MeV, where a significant amount 
of the 2p2h component is predicted [8,9]. 

In 1997, Wakasa et al. [10] accurately measured the ^°Zr(p, n) spectra at 295 
MeV, the energy at which spin-fiip cross sections are large, distortion effects 
are minimal [11] and therefore the characteristic shapes of the angular dis- 
tributions for each angular momentum transfer (AL) are most distinct. They 
successfully identified the GT strength in the continuum region through mul- 
tipole decomposition (MD) analysis which extracted the AL = component 
from the cross sections [10]. They obtained a GT quenching factor, defined as 

Q = --y— of 0.90 ± 0.05, where the error is due to the uncertainty of 

the MD analysis [10]. As discussed in Ref. [10], the main source of the system- 
atic uncertainties is the overall normalization, i.e., the GT unit cross section 
and the S13+ value. The uncertainty of the GT unit cross section amounts to 
16%, but will be reduced to ~ 5% by the ongoing systematic analysis of new 
(Tgt data at 295 MeV [12]. On the other hand, the uncertainty in the Sp+ 
value is difficult to properly assess. A S13+ value of 1.0 ± 0.3 used in Ref. [10] 
was obtained from a similar, but simpler MD analysis of ^°Zr(n,p) data at 
198 MeV by Raywood et al. [13]. Raywood et al. found a significant amount 
of the monopole (AL = 0) cross sections in the continuum in the region of 
Ex — S 31 MeV, which corresponds to a GT strength of ~5, but attributed 
all to the isovector spin-monopole (IVSM) strength [13]. Thus, the quenching 
factor was subject to uncertainties of both the f3~^ strength in the continuum 
and contributions of the IVSM component. To reduce those systematic uncer- 
tainties, it is essential to have accurate {n,p) data at the same energy as the 
{p, n) data, and to perform consistent analyses on both sets of data. In this 
letter the measurement of the ^'^Zr(n,p)^°Y reaction at 293 MeV is reported. 
Performing a consistent analysis on both the (p, n) and (n, p) data, we have 
derived a reliable GT quenching factor. 

The measurement was performed with the (n,p) facility [14] at the Re- 
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search Center for Nuclear Physics (RCNP). A schematic view of the {n,p) 
facihty is shown in Fig. 1. A nearly mono-energetic neutron beam was pro- 
duced by the ''Li(p, n) reaction at 295 MeV. The primary proton beam, after 
going through the ^Li target, was deflected away by the clearing magnet to a 
Faraday cup in the floor. The typical beam intensity was 450 nA and the thick- 
ness of the '''Li target was 320 mg/cm^. About 2 x 10^/s neutrons bombarded 
the target area of 30^^ x 20"^ mm^ located 95 cm downstream from the ''Li 
target. Three ^"^Zr targets with thicknesses of 485, 233, and 215 mg/cm^ and 
a polyethylene (CH2) target with a thickness of 46 mg/cm^ were mounted in 
a multiwire drift chamber (MWDC). Wire planes placed between the targets 
detected outgoing protons and enabled us to determine the target in which 
the reaction had occurred. Charged particles coming from the beam line were 
rejected by a veto scintillator with a thickness of 1 mm. The scattering angle of 
the (n, p) reaction was determined by the information from the target MWDC 
and another MWDC installed at the entrance of the Large Acceptance Spec- 
trometer (LAS) [15]. The outgoing protons were momentum-analyzed by LAS 
and were detected by the focal plane detectors [16]. The number of ^H(n,p) 
events from the CH2 target was compared to the SAID [17] calculated ^H(n,p) 
cross sections for normalization of the neutron beam flux. Blank target data 
were also taken for background subtraction. 

We have obtained double differential cross sections up to 70 MeV excitation 
energy over an angular range of 0°-12° in the laboratory frame. The data have 
been analyzed in 1-dcgrce bins. The ^'^Zr(n,p) spectra at three of the twelve 
angles are shown in the left panel in Fig. 2 by the solid dots. The overall 
energy resolution is 1.5 MeV, mainly originating from the target thicknesses 
and the energy spread of the beam. The angular resolution is 10 mr which 
is dominated by multiple scattering effects in the ^''Zr targets. In addition to 
the statistical uncertainty of ~ 2% per 2-MeV excitation energy bin, there 
is a systematic uncertainty of 5%, where the main contributions are that of 
target thicknesses (4%) and the angular distribution of the n+p cross section 
taken from the phase-shift analysis (2%). The right panel in Fig. 2 shows 
the ^'^Zr(p, n) spectra [10]. At 0° the {p,n) cross sections arc larger than the 
(n, p) cross sections not only in the GT resonance region but also in the high 
excitation energy region of = 70 MeV due to excess neutrons in ^'^Zr. 

The MD analyses were performed on the (p, n) [10] and (n,p) excitation 
energy spectra to obtain GT strengths. Details of the MD analysis are given 
in Ref. [10]. For each excitation energy bin between MeV and 70 MeV, 
the experimentally obtained angular distribution a^^^{9cin, Er^) has been fit- 
ted using the least-squares method with the linear combination of calculated 
distributions, 

Cr'^^^'^iOcm: E^) = ^ ajTvapf^^n (^cm, E^) , (1) 

where the variables ajn are fitting coefficients all of which have positive val- 
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ues. The calculated angular distributions for each spin and parity transfer 
^ph^j^i^cm, E,j) have been obtained using the distorted wave impulse approxi- 
mation (DWIA) calculations described below. 

The DWIA calculations were performed with the computer code DW81 [18] 
for the following ^ transfers: 1+(AL = 0), 0", 1", 2-(AL = 1), 3+(AL = 2), 
and 4~(AL = 3). The one-body transition densities were calculated from pure 
Iplh configurations. The (1(77/2, l(/(7/2) ^'^^ (15^9/2, l5'9/2) configurations were 
used to calculate the GT transitions in the analyses of both the (p, n) and 
the (n,p) spectra. For the transitions with AL > 1 in the (p, n) channel, the 
active proton particles were restricted to the l5'9/2, l5'7/2, 2^5/2, 2^3/2, l/iii/2, 
or 3s 1/2 shells, while the active neutron holes were restricted to the 15^9/2, 2pi/2, 
2p3/2, 1/5/2, or 1/7/2 shells by assuming "^^Ca to be the core. In the analysis 
of the spectra, the active neutron particles were restricted to the lg-j/2i 

2ci5/2, 2(^3/2, l/iii/2, or 3si/2 shells, while the active proton holes to 2pi/2, 2p3/2, 
1/5/2, or 1/7/2. The optical model potential (OMP) parameters for proton were 
taken from Ref. [19]. The OMP parameters for neutron were also taken from 
Ref. [19], but without the Coulomb term. The effective A^A^ interaction was 
taken from the t-matrix parameterization of the free A'^A'^ interaction by Franey 
and Love at 325 MeV [1]. It should be noted that DWIA calculations using 
this parameter set better reproduce the polarization transfer L'jv7v(0°) for the 
^°Zr(p, n) reaction than those at 270 MeV [10]. The radial wave functions 
were generated from a Woods-Saxon (WS) potential [20], adjusting the depth 
of central potential Vo to reproduce the binding energies [21,22,23,24]. The 
unbound particle states were assumed to have a shallow binding energy to 
simplify the calculations. For a given J'^ transfer, the shapes of the angular 
distributions depend on the Iplh configurations. Thus, all combinations of 
Iplh configurations were examined by the ^^-minimization program and the 
optimal combination was obtained for each excitation energy bin. 

Results of the MD analyses are shown in Fig. 2. The obtained AL = 
component in the (p, n) spectra has a large contribution not only in the GTGR 
region, but also in the high excitation energy region up to 50 MeV [10]. The 
AL = component of the cross section, (7al=o(?, <^), is proportional to the 
GT strength S(GT) [5] such that, 

aAL=o{q.uj) = crGTF{q,uj)B{GT), (2) 

where (Tgt is the GT unit cross section [10] and F{q,u!) is the kinematical 
correction factor [25] . The GT unit cross section has been determined so that 
the Sp- value up to the GTGR region of < 16 MeV becomes 18.3 ± 3.0 [5] 
as described in Ref. [10]. The obtained value is o-qt = 3.5 ± 0.6 mb/sr, which 
is consistent with 3.6 ± 0.6 mb/sr, the value used in Ref. [10]. 

The strength distributions are shown in Fig. 3. Here the contribution from 
the isobaric analogue state (IAS) at 5.1 MeV, corresponding to 0.7 ±0.1 in GT 
unit [10], is already subtracted. The strength is denoted as S(GT-I-IVSM) 
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because it contains the IVSM component [13]. The error bars are the ilcr 
confidence hmits obtained by a Monte-Carlo simulation, where the mini- 
mization is performed for synthetic data sets generated by replacing the actual 
data set in accordance with the statistical errors [26] . The MD analysis of the 
(p, n) spectra becomes unstable above 50 MeV excitation [10] and reanalysis 
with larger bin width has not improved the stability. Therefore we set the 
upper limit of the excitation energy to 50 MeV. The integrated strength thus 
obtained is 'S'(p_„);gt+ivsm = 33.5 ± 0.6(stat.) ± 0.4(MD) ± 4.7((7gt) up to 
50 MeV excitation energy, neglecting the uncertainty in the subtraction of the 
IAS contribution. The first error is the statistical error of the MD analysis. The 
systematic uncertainty due to the input parameters for DWIA calculations has 
been evaluated by using the wave functions generated from a potential by the 
relativistic Hartree approach [27] or with other OMPs [11,27,28], and it is 
estimated to be ±0.4. The third error reflects the error in the GT unit cross 
section. The distribution of strength in Fig. 3 is shifted by +18 MeV, 
accounting for the Coulomb displacement energy and the nuclear mass differ- 
ence. The strength integrated up to 32 MeV excitation of ^"^Y, or = 50 MeV 
in Fig. 3, is %,p);GT+ivsM = 5.4 ± 0.4(stat.) ± 0.3(MD) ± 0.9(aGT). 

The curves in Fig. 3 show the theoretical predictions of GT strength 
distribution by employing the dressed particle random phase approximation 
(DRPA) model [29] , folded by a Gaussian distribution to simulate the energy 
resolution of the measurement. The agreement between the experiment and 
the theory is excellent, except in the excitation energy region around 30-40 
MeV, where the IVSM resonance is important [13,30]. The IVSM resonance is 
the 2hu! excitation via the r^aT± operator and has the same spin and parity 
transfer as the GT resonance. Since the angular distribution of the IVSM tran- 
sition has a forward peaking shape similar to that of the GT transition [31,32], 
the present MD analysis cannot discriminate these two components. 

Hamamoto and Sagawa studied the IVSM modes for both the ^°Zr(p, n) 
and ^°Zr(n,p) reactions [30]. They pointed out that the response function to 
the operator r^(Tr± calculated for the 2?iuj excitation contains a small GT 
component because of the difference of the neutron and proton one-particle 
wave functions with the same quantum numbers. This GT component has to 
be subtracted to obtain the pure IVSM strength [30]. Although interference 
exists between the GT and the IVSM modes, the cross section associated with 
the pure IVSM component is estimated in this work and its contribution is 
subtracted incoherently from each spectrum since the distribution of the GT 
strength in the IVSM resonance region is unknown. The DWIA calculations 
have been performed by assuming that the strengths are fully exhausted in the 
state with central energies reported to be 35 MeV in the (j),n) spectrum [32], 
corresponding to 19 MeV excitation of ^°Y, and — 37 MeV in Fig. 3, in the 
(n,p) spectrum [30]. The transition densities are obtained by the procedure 
of Conde et al. [33]. The small GT components are explicitly eliminated by 
modifying the radial wave functions of the final states. The obtained IVSM 
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strengths are within 10% of the theoretical prediction [30]. The calculated 
IVSM cross sections at 0° for the {p,n) and the {n,p) channels are 6.9 ± 
1.5 mb/sr (4.2 ± 0.9 GT units with = 0.47) and 5.3 ± 0.6 mb/sr 

(2.5 ± 0.3 GT units with F{q,uj) = 0.61), respectively. The uncertainties are 
mainly due to the choice of OMP parameters. By subtracting these values from 
^GT+ivsM, the actual GT strengths of Sg- = 29.3 ± 0.5(stat.) ± 0.4(MD) ± 
0.9(IVSM) ± 4.7(aGT) and Sp+ = 2.9 ± 0.4(stat.) ± 0.3(MD) ± 0.3(IVSM) ± 
0.5((Tgt) have been obtained. This new S/s- value is consistent with but slightly 
higher than the one reported in Ref. [10]. The Sp+ value agrees well with the 
DRPA prediction of Sf)+ = 3.2 [29]. 

The present strength may be compared with previously reported re- 
sults that employed similar MD techniques on the ^^ZT{n,p) spectra. Conde et 
al. [33] obtained Sfi+ = 1.7±0.2 up to 10 MeV excitation at T„ = 98 MeV while 
Raywood et al [13] obtained 1.0 ± 0.3 up to - 8 MeV at T„ = 198 MeV. The 
errors given above do not include the uncertainty in the GT unit cross section. 
The /3+ strengths obtained in this work up to 8 MeV and 10 MeV excitation are 
Sp+ = 0.4±0.1(stat.)±0.1(CTGT) and 0.7 ± 0.1(stat.) ± O.l(aGT), respectively. 
These 3^+ values are significantly smaller than those obtained by Raywood et 
al. or by Conde et al. We note that the analysis of the T„ = 98 MeV data may 
suffer from ambiguities in the reaction mechanism due to multi-step processes 
while the MD analysis at 200 MeV may suffer from the ambiguity due to the 
distortion effects which are larger than those at 300 MeV. 

The consistent analyses of both (p.n) and (n.p) spectra yield a quenching 
factor of Q = 0.88±0.02(stat.) ±0.05(syst.) ±0.01(MD) ±0.02(IVSM), where 
the systematic uncertainty of the normalization in the cross section (5%) is 
also indicated. An uncertainty in the GT unit cross section of 16% is not 
included. It should be noted that since the errors are correlated, the combined 
systematic errors are smaller than their geometric sum. Thus quenching of GT 
strength up to 50 MeV due to coupling between ph and Ah states becomes 
significantly smaller when 2p2h contributions are properly accounted for [9,10]. 

The interpretation of this small quenching in terms of the short range 
correlation in nuclei is particularly interesting. Assuming that the missing GT 
strength of ~ 10% is attributed to the A excitation, one can derive the Landau- 
Migdal (LM) parameter g^^ which describes the short range correlations for 
ph Ah transitions in the vr + p + model [34]. The g^^ value deduced by 
using an RPA model in ph and Ah spaces [34] with the Chew-Low [35] coupling 
constant (/a/Ztt = 2) is shown in Fig. 4 as a function of Q. If we take here 
Q — 0.88 ± 0.06, combining the uncertainty of MD analysis and that of IVSM 
contribution in quadrature, then we derive g^^ = 0.18 ±0.09 assuming g'^^ ~ 
0.5 [36]. Arima et al. have examined the finite size effects of the ^"Zr nucleus by 
taking the finite range interaction due to tt- or p-exchange into account [37] . If 
their argument is employed here, the ^^^a value increases by 0.07 for the same 
Q. Therefore it is reasonable to assume the g'^^ value to be 0.25 ± 0.09, which 
is significantly smaller than that assuming the universality relation of the LM 
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parameters, i.e., g^^ = g^^ = g'^^ = 0.6-0.8 [38]. This result is consistent 
with those obtained using the coupled channel G- matrix calculations [36,39]. 
An important consequence of such a small ^i^a value is the enhancement of 
the pion correlation in nuclei. According to the RPA prediction by Tatsumi 
et al. [40] the critical density of the pion condensation becomes (1.9 ± 0.3)po, 
just a half of that predicted by employing the universality assumption. 

In summary, we measured the double differential cross sections between 
0°-12° for the ^°Zr(n,p) reaction at 293 MeV in an excitation energy region 
of 0-70 McV to study the GT strengths in the continuum. The MD tech- 
nique was applied to the measured cross sections to extract the AL = cross 
section in the continuum. After subtracting the IVSM contribution a total GT 

strength of ^ = 2.9±0.4(stat.)±0.3(syst.)±0.3(MD)±0.3(IVSM)±0.5(aGT) 
up to 32 MeV excitation was obtained. A revised and consistent analysis of 
the (p, n) [10] and {n,p) reaction data from ^°Zr yield a reliable quenching 
factor of Q = 0.88 ± 0.06, not including the uncertainty of the GT unit cross 
section of 16%. The error includes the uncertainty of the estimated IVSM 
contribiition as well as error in the data and the MD analysis. This work is 
the first attempt to deduce the GT quenching factor by accurately taking into 
account GT strengths in the continuum region. 
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Fig. 1. A schematic view of the {n,p) facihty at RCNP. 




(MeV) 

♦ data ^MAZ,=?0 :raAL = l EZiAL = 2 EsaAZ,=3 



Fig. 2. Double differential cross sections for the ^''Zr(n,p) (left panel) [10] and 
^°Zr(p, n) (right panel) reactions. The histograms show the results of the MD anal- 
yses. 
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Fig. 3. GT plus IVSM strength distributions obtained by the MD analysis of the 

^^ZT{p,n) and ^^ZT{n,p) reactions (in GT unit). The ^°Zr(n,p) spectrum is shifted 
by +18 MeV. The curves are taken from Ref. [29]. The energy regions of IVSM 
excitation are indicated by braces. See text for details. 
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Fig. 4. LM parameter g'^^ in Chew-Low model estimated as a function of quenching 
factor Q. The finite size effect of ^^Zi nucleus is not taken into account. 
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